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Abstract 

Diffraction of a surface wave on a rectangular wedge with impedance faces is studied using the 
Sommerfeld-Malyuzhinets technique. An analog of Landau's bypass rule in the theory of plasma waves 
is introduced for selection of a correct branch of the Sommerfeld integral, and the exact solution is given in 
terms of imaginary error function. The formula derived is valid both in the near-field and far-wave zones. 
It is shown that a diffracted surface wave is completely scattered into freely propagating electromagnetic 
waves and neither reflected nor transmitted surface waves are generated in case of bare metals which have 
positive real part of surface impedance. The scattered waves propagate predominantly at a grazing angle 
along the direction of propagation of the incident surface wave and mainly in the upper hemisphere regard- 
ing the wedge face. The profile of radiated intensity is nonmonotonic and does not resemble the surface 
wave profile which exponentially evanesces with the distance from the wedge face. 

Comparison with experiments carried out in the terahertz spectral range at Novosibirsk free electron laser 
has shown a good agreement of the theory and the experiments. 

PACS numbers: 42.25.Bs; 42.25.Fx; 42.25.Gy 
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I. INTRODUCTION 



Plasmonics is now a very rapidly developing field of activity, and advancement to the mid- and 
far-infrared (terahertz) ranges of frequencies is one of the main streams in photonics [1, 2]. Our 
interest in this subject has arisen from discussions of the results of experiments on propagation 
of Surface Plasmon Polaritons (SPPs) along gold-ZnS-air interfaces [3], which were carried out 
at Novosibirsk Free Electron Laser facility [4] in the terahertz spectral range. For SPPs not to 
be disturbed by any material probe, plasmon-polariton characteristics were studied indirectly via 
sensing of the electromagnetic field in the space behind the tail facet of samples. In the exper- 
iments, the profile of radiated intensity is non-monotonic and does not resemble a surface wave 
profile which exponentially evanesces with the distance from the wedge face. A maximum of the 
radiation intensity is observed at some distance from the surface plane whereas a maximum of 
the surface wave field should be located at the wedge surface. These discrepancies had initiated 
search for theoretical explanation which has led us to the theory developed by G. D. Malyuzhinets 
in the 1950s. 

The Sommerfeld-Malyuzhinets theory is known to provide a powerful method for exact so- 
lution to specific optical problems [5-7]. In particular, G.D. Malyuzhinets solved the problem 
of diffraction of a surface wave by an impedance wedge, considering it as a special case of a 
plane wave propagation at the Brewster angle [8]. He derived some general relations but he did 
not investigate the properties of integrand functions near the saddle points which is necessary for 
computation of the integrals. A few books and review papers [5-7] of varying elaboration provide 
a systematic introduction to the theory by G.D. Malyuzhinets. Some earlier approaches to the 
problem of diffraction by wedges and screens are reviewed in [9, 10]. We would like to add a few 
references to this list [11-21]. 

V. Zon in Ref. [9, 10] applied the Sommerfeld-Malyuzhinets technique to the diffraction of 
a surface wave on conducting wedge with small surface impedance but she did not succeed in 
obtaining the final solution. We do that in this paper. We also conclude that her calculations of 
the reflection and transmission coefficients for the surface wave are incorrect. To obtain correct 
results we introduce a sort of Landau's bypass rule known in the theory of plasma waves [22]. 

In Section II, we briefly recall the properties of surface waves to a minimal extent necessary 
for understanding of the subsequent calculations and justify the use of the Leontovich boundary 
conditions at the wedge faces; in particular, we give (without derivation) a formula for the sur- 
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face impedance of a metal substrate coated with a thin dielectric film. The foundations of the 
Sommerfeld-Malyuzhinets technique are expounded in Section III. The diffracted field is calcu- 
lated in Sections IV, V, and VI. A brief description of the experimental techniques and compar- 
ison of experimental results with theoretical predictions are given in Section VII. Section VIII 
concludes the paper and summarizes our results. 



H. SURFACE PLASMONS ON A PLANE 

Let a metal with the complex permittivity s occupy the lower half-space y < 0. Since a TE 
wave (the electric field is transverse to the plane of incidence) has no surface branch, we consider 
a TM wave in the upper half-space y > 0, localized near the metal-air interface and propagating in 
the x-direction: 

E(r, t) = E () exp [-xy + i(kx - cot)} , E = {E Qx , E Qy , 0}, 
B(r, t) = B exp [-xy + i(kx - cot)] , B = {0, 0, B 0z ). 

It is characterized by the circular frequency co, wavenumber k, and attenuation constant x with a 
positive real part, Rex > 0. Similar fields (with a negative attenuation constant, Rex' < 0) could 
be written for the lower half-space but they are needed only for derivation of the dispersion law 
for the surface wave: 

k = k S^. (2) 

V 1 + £ 

Here 

k = co/c, x 2 = k 2 - 

the Gaussian system of units is used, and the derivation can be found elsewhere (see [23, 24]). 
Eq. (1) can be interpreted as an evanescent plane wave (also known as the Zenneck wave [25, 26]) 
propagating at a complex angles such that 

k x = ^ocos^ = k, k y = -ko sin^ = ix, x = iko sin^. (3) 

Using dispersion relation of Eq. (2) yields 



In a fictitious case of real e, the surface wave exists provided that e < -1. For complex values of e 
there is always a solution with Im^ < which exponentially decreases with the distance from the 



boundary of the metal. However, speaking about surface waves it makes sense only if the damping 
length \/lmk x is sufficiently long, i.e. \lmk x \ «: |Refc x |. This condition is certainly satisfied if 
\s\ » l.Then k x ~ k (1 - 1 /2s) and the imaginary part of k x is automatically small no matter how 
large or small the imaginary part of s is. 

Keeping the terahertz radiation at Novosibirsk Free Electron Laser [4] in mind, we will take as 
an example the value of permittivity 

s = -103260 + / 310810 

for gold at a frequency corresponding to the wavelength A = 140 microns [27]. Then 

X = 0.001022 - / 0.001417, 
k x = (44880 + / 0.065) m" 1 , 
x = (64 + /46)m _1 . 

The components of the electric and magnetic field at the metal surface are related through the 

boundary condition, which can be derived from the Maxwell equation 

„ IdE 
VxB = -— . 

c dt 

Its ^-component 

-kB 0z = -ik E 0x 

together with Eq. (3) yields the Leontovich boundary condition E 0x = sin^5 0z -In the general case, 
the Leontovich boundary condition relates the tangential component of the electric field E T with 
that of the magnetic field. It is usually written in the form E T = £ [n x B] , where the unit vector 
n is directed along the outward normal to the surface of the metal, and the parameter £ is called 
(dimensionless) surface impedance. 

In Leontovich's theory, £ is assumed to be a function of frequency, which depends only on the 
conductor material but neither on the incident angle nor the type of incident wave. Using Fresnel's 
formulae one can readily check that 




for a plane TM wave impinging against the metal surface under the angle 9 to the normal. A 
propagating TE wave is characterized by a different surface impedance 



and 



1 



(5) 



£ = sin* = 



l+e 



for the surface wave (1). In any case, 



£- 1/V£ 



provided that inequation \s\ » 1 holds. We see that for large \s\ the surface impedance £ is approx- 
imately independent on parameters of the waves, which justifies using the Leontovich boundary 
condition on a metal surface instead of looking for complete solution for wave propagation inside 
a metal. 

Below we will write the Leontovich boundary condition in the following form 



assuming that the formal parameter* is related with the surface impedance by Eq. (5). As shown 
above,* can be considered as a complex angle of propagation of surface wave. For the above cited 
parameters of golden surface, 



Coating of metallic surfaces strongly affects the magnitude of £. For a metal coated with a thin 
film with a small width d <sc 2n/k and permittivity » 1, we derived the following expression 
for the surface impedance: 



Properties of the film becomes dominating for d > 1 /k yfs. For example, ZnS film with d = 
0.75//m and 

s d = 8.7 + i 0.059 

over a gold substrate has the impedance 



E T = sin* [n x B] , 



(6) 



£ = 0.00102-i 0.00142. 



(7) 




(8) 



£^0.00106-/0.0312 



(9) 



with an imaginary part more than 20 times larger than that of pure gold. 



III. SOMMERFELD-MALYUZHINETS THEORY 



To begin with, let consider a plane wave of unit amplitude incident on a wedge at the angle 
6 = 8 > 0, as shown in Fig. 1. In polar coordinates r and 6 perpendicular to the edge of the 
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Figure 1. (color online) The geometry of the problem in cylindrical coordinates. The angles 9 and <1> are 
measured from the bisector of the wedge. 

wedge, the wave is written as 

5; nc = exp [-/fc o rcos(0 - 6q)] . (10) 

Let the wedge occupy the region \6\ > O. Within the framework of geometrical optics, the upper 
face of the wedge reflects the incident wave at the angle 6 = 2<D - 6 - n. One more n should be 
subtracted from this value for taking the reversal of the direction of propagation into account, thus 
the field of the reflected wave is derived from Eq. (10) with the substitution 6 — » 2$ - 6 - In — » 
20 - 6 : 

5 refl oc exp [-ik r cos(6> + 9 Q - 20)] . (11) 

Below we restrict ourselves to the case of TM waves with the transversal component of mag- 
netic field B z = B. The reason for such choice is that a surface wave may have no TE polarization. 
In the vacuum region, |0| < O, the TM wave obeys the Helmholtz equation 

V 2 B + k 2 Q B = 0. (12) 

The non-zero components of the electric field in the TM wave, 



are expressed through the derivatives of B. The Leontovich boundary conditions in Eq. (6) with 
the surface impedance of Eq. (5) on the conducting faces of the wedge at 6 = +<D now take the 
following form 

IdB 

-— = +ik smx±B, (14) 
r 06 

where the + and - signs label the quantities related to the faces of the wedge located at the angle 
= +0 and 6 = -O, respectively. 

Heuristic considerations enable construction of a solution to Eq. (12) with given boundary 
conditions [7]. According to Sommerfeld [28], the solution is sought in the form of a superposition 
of plane waves 

l r 

B{r,9) = — exp (-ik r cos p) s{6 - p) dp, (15) 
2ni J 

y 

where the integration is carried out over the contour y in the complex p-plane, k r > and the 
kernel s(6 - p) = s(9 - p, 6 , <D) is analytic outside the real axis. The shadowed areas in Figure 2 
indicate the regions of the complex p-plane where Imcos p < and the exponent of Eq. (15) tends 
to zero as kr —> <x>. 

Following Sommerfeld, we choose y in the form of two loops. Then, the integral in Eq. (15) is 
a solution to Helmholtz's equation, provided that the two loops, y + and y_, are shifted upwards and 
downwards correspondingly by the distance |lm(0-p)| > |Im(y)|; thus any singularities (see below) 
of the function s{6 - p) are located between these loops within the strip -|Im(^)| < Im(# - p) < 
\lm{x)\- Choosing the integration contour in the form of two symmetrical loops guaranties that the 
asymptote of the integral in Eq. (15) at k r —> oo contains no converging cylindrical waves except 
for the incident wave of Eq. (10), and the incident wave is introduced into s through given a pole 
at p = 8o as explained in [7]. 

Let now write functional equations for s. For the Leontovich boundary conditions of Eq. (12), 
we have 



exp (-ikorcos, p) 



- 5'(±0 - p) + ik sinx± s(±0 - p) 
r 



dp = 0, (16) 



where s' = ds/dp = -ds/d9. Integrating the first term by parts yields the following equations 

exp (-ik r cos p) [- sin p + sin^ ± ] s(±<£> - p) dp = 0. (17) 



r 

Since the contour y = y + U y_ consists of two loops situated symmetrically to the point p = 0, 
these equations are true if the kernel s is an even function of p. This leads us to the Malyuzhinets 




Figure 2. (color online) Plane of the complex angle p with the Sommerfeld double loops y + U y_ (solid 
line) and the steepest descent paths go U gi U g2 (dashed line). 

functional equations [8]: 

[sin p + sin^+] ^(p + <D) = [- sin p + sin^ + ] s(-p + O), 
[sin p - sin^„] s(p - <D) = [- sin p - sin^_] s(-p - O). 
With rather non-trivial calculations one can verify [7] that the function 



(18) 



(19) 



where 



*(v,Ar) = ^® (v + ® + ±tt ^ (v + O - Itt +x) , 

1 r°° cosh(w?7) - 1 
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(w) = exp 



d?7 



(20) 
(21) 



7 cosh ^7T77 ) sinh(2®77) 

supplies a particular solution for s(u) to the functional equations (18). The Malyuzhinets function 
iff® is regular in the strip |Re(p)| < \it + 2®, where the integral in Eq. (21) is even and satisfies the 
following functional relation 



iff® (w + 2<D) /iff® (w - 2<D) = cot (jw + \n) . 
With Eq. (22) true this function extends beyond the indicated band. 



(22) 



For O = mn/An, where m/n is an irreducible rational number, the derivative din fa/dp of In 
can be written as a finite sum of trigonometric functions [8]. In particular, 

. 4 (w - n\ (w + n\ (w\ 4 (w\ 1 (w\ 

(w) = 3 cos (— j cos (_) sec (- j = - cos (-) - - sec (-) (23) 

for O = |tt, which corresponds to a right angle wedge. 
Looking for a general solution in the form 

s(u) = cr^^oW/Wo), (24) 
we deduce from Eqs. (18) that the function cr(u) obeys the following equations 

cr(n + O) = cr(-p + 0), 

(25) 

<j{p - 0) = cr(-p - 0). 

For description of an incident wave with unit amplitude, cr(6 - p) should have a simple pole at 
p = 9q with unit residue. By virtue of Eqs. (25), <x is symmetric about points p = +0, and 
therefore its system of poles should satisfy the same symmetry relationships. If the pole p = 9 
is successively reflected about points p = and p = -0, one obtains for s(Q - p) a lattice of 
poles with the period 40 at p = 6 + 4n0 (n = 0, +l,+2, etc.) and a similar grid at points 
p = 20 - # + 4n0. It is possible to guess a function with these properties. It was found by 
G.D. Malyuzhinets [29]: 



jucos0u# ) 

cr(u) ~ 



where p = 7r/20. For = \ n 



sinQuu) - sin(//#o) ' 
o-( u ) = , v ' o , . (26) 



sin(f)-sin(f) 



IV. APPROXIMATION OF GEOMETRICAL OPTICS 

At a first glance, the Sommerfeld-Malyuzhinets theory might be sought of as a purely math- 
ematical trick. However, it indeed describes a real physical phenomenon. To prove that, in this 
Section we show how an approximation of geometrical optics can be retrieved from this theory 
which provides an exact solution of the problem of vector field diffraction, although for a very 
limited set of problems related to the diffraction of electromagnetic waves on a conducting wedge. 
To avoid possible misunderstanding regarding the role of geometrical optics in the theory to be 
presented in the next three sections, it is worthy to emphasize that the field which will be derived 



in the current section is merely a part of a complete solution. From mathematical point of view, 
it represents a contribution of the poles of the integrand core s(9 - p) in Eq. (15) whereas a com- 
plete solution (49) includes also a contribution of the saddle points p = +n. The separation of the 
contribution of the poles and saddle points can be approximately made for the distances r from the 
wedge top larger than the wavelength, k r » 1. 

Following G.D. Malyuzhinets [8], we replace the contour y = y + U y_ with the 3-loop contour 
go U gi U g 2 shown in Fig. 2 with the dashed line, in order to compute the Sommerfeld integral in 
Eq. (15): 

B(r, ff) = J- f exp (-ik r cos p) s{9 - p) dp. (27) 
2m J 

The contours go, gi and g 2 are obtained by deformation of the contours y + and y_ into the steepest 
descent paths through the saddle points p = -n and p = n of the exponent exp(-/fc rcos p) in the 
integrand. The new contour go is formed from the middle parts of the contours y + and y_, whereas 
the contours gi and g 2 are composed by remaining end parts of y + and y_. As will be shown in 
Section V, the contribution of the saddle point describes the diffracted field at large distances (in 
the wave zone). 

The integral over the contour g is reduced to the sum of the residues in the poles of the inte- 
grand encircled by g within the strip 

-n < Re(p) < n. (28) 

These are the poles to be crossed during the deformation of the original contour of integration 
yi U y 2 to the steepest descent path contour g\ U g 2 . The number of such poles depends on the 
values of <D, 6, and 8 . 

Let consider first the poles of the function cr in Eq. (24). There are an infinite number of such 
poles located at the points 

p n = 9-(-l) n e -2n^, n = 0,±l,.... (29) 

However only three of them, namely 

p = 9 - 9o, Pl = 9 + 9 - 20, p_ l = 9 + 9 + 20, (30) 

have chance to fall into the region (28) if O > \n. 

The residue at the pole p stands for the wave exp[-zfcor cos(#-# )] , which is readily recognized 
as the incident wave. It enters into the result of integration if -n < 9 - 9 < n. Values of 9 outside 
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the interval -n + 80 < < n + 9 and within the interval -0 < < represent the regions of 
geometrical shadow for the incident wave. 

The pole p x gives rise to the wave exp [ -z'fc r cos (# + #0 - 2<D)] = exp[z'/c r cos(# + 6 - 20 + n)], 
reflected from the upper face 6 = of the wedge (see Fig. 1). It enters into the result of integration 
if -n < 6 + 0„ - 20 < n, i.e. 20 - 6 - 6 Q < n. 

Finally, the pole p_i describes the wave exp[-z'/c rcos(# + 9 + 20)], reflected from the lower 
face 6 = -O of the wedge. It enters into the result if 20 + 6 + 6 < n. 

Thus, the geometrical optics yields a solution constituted by 'pieces' of plane waves. At the 
ends of these pieces represented by light-shadow boundaries the solution vanishes jumpwise to 
zero in the shadow region, which means that the geometric-optical solution is discontinuous. These 
discontinuities of the optical fields will be eliminated with the addition of the diffracted fields in 
Section V. 

The surface wave, grazing towards the wedge top along the upper face, as shown in Fig. 1, has 
a complex propagation angle: 

0o = 0-* + . (3D 

It cannot be reflected from the lower face of the wedge if > \n as can be readily deduced from 
the above treatment, and only two poles, namely 

p = 6-®+x + , Pl =9-0>-x + , (32) 

can fall into the region given by inequation (28). This occurs if 

- 7T + Re^ + < 6 < 0, (33) 

where the upper and the lower signs stand for p and pi, correspondingly. 

Due to the factor l/*P o (0o) in Eq. (24), the residue about the pole p is evaluated as a wave of 
unit amplitude: 

B (r, 6) = exp [-ik r cos(6> - + x + )] ■ (34) 

It describes the incident surface wave propagating along the upper face to the top edge of the 
wedge. According to the inequation (33), this wave does not penetrate into the shadow region 
6 < - n if Re^f + is sufficiently small. We will see in Section V that the transitional region near 
the formal boundary 9 = - n - Re^ + of the shadow region is described by a simple function, 
which includes the contributions from both the pole p and the saddle point p = -n. 
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The contribution of the pole p x is exactly zero, which can be seen from the fact that formally 
calculated residue ¥0(8 - Pi) contains the multiplier 1/^(20 + n/2), which is zero. Since the 
contributions of both the p x and p_\ poles are zero, the wedge does reflect surface waves (this 
assertion will be clarified in Section VI). 

Let us now proceed to the poles of the function ¥0(6 - p). To simplify our task, below we 
restrict ourselves to the case of a right angle wedge with <D = jn. Then, *F O (0 - p) is expressed 
through trigonometric functions with the aid of Eq. (23), and the poles can be found from the 
following equation 



cos 







cos 









X cos 



5n 
T 



1 L 57r 
z [°-x+-p + T 



cos 



= 0. (35) 



They obey the conditions of Eq. (28) if 



(36) 



<6<-\n -Rex- 



Vox bare metals, Re^ ± is positive in the optical and infrared ranges of frequencies. Hence, the 
inequations (36) cannot be satisfied, and the function ¥ cannot have poles inside the contour go- 
However the zero point 

p + = 9-x + -l* (37) 

of the factor cos [| (9 ~x+ - P + x)] approaches the first saddle point p = —n at 6 = O = \n. 
Similarly, the zero point 

p-=e+x- + }n (38) 

of the multiplier cos ^ (d +x~ ~ P ~ lfj\ is located near the second saddle point p = +n at 6 = 
-\n. Although one might expect that these poles strongly affect the diffracted field near the wedge 
faces at 6 = calculations in Section VI do not confirm these fears. 



V. DIFFRACTED FIELDS 

Evaluating the integral in Eq. (27) along the contours g\ and g 2 in Fig. 2 yields radiated (freely 
propagating) electromagnetic fields. Those contours pass through the saddle points p = —n and 
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p = n, respectively. To use the saddle-point method of integration, we expand the exponent in the 
integrand of Eq. (27) about the saddle points so that 

1 -i 

p = +7T + —— t. 

^/k^r 

Putting the expansion 

expi-ik^r cos p) exp[ik r - t 2 ] 

in Eq. (27) reveals that, due to the factor exp(-? 2 ), the main contribution to the integral over the 
contours gi and g 2 comes from the neighborhood of the saddle points. Turning to the integration 
in the variable t and taking into consideration the direction of integration over the contours gi and 
g2, we transform the respective integrals to the following form: 

CO 

Bl ( r ,d) = J-. f s(9-p)e- ikorcosp dp = -^-. f s\e + n-^^t\e ikQr - t2 ^-dt, (39a) 
2m J ■ 2m J \ ^Jk~? j dt 

oo 

B 2 (r, ff) = — f s(9- p) Q- ik ° r ™Pdp = +— f sle-n- -^4 t) — dt. (39b) 

2m J " 2m J \ yp^r J dt 

g2 -00 

In the wave zone that corresponds to the limit k r —> oo, the pre-exponential factor s(9 - p) can be 
substituted with s(9 + n), which yields 

gilcQr+in/4 

B l (r,9) = —=s(9 + n), (40) 
y2nkor 

for the contour gi and 

B 2 (r,9) = -—= S (9-n) (41) 
y2nk r 

for the contour g 2 . 

Since the pole p = 9 - |tt +x+ approaches the saddle point p = —n when 9 tends to the 
function s(9 + n) has a narrow peak near 9 = -\n, which is much bigger than s{9 - n) as shown in 
Fig. 3. Hence, the field B 2 can be neglected near the peak, where the amplitude of the diffracted 
field is 

i exp (ik Q r + jin) 

By - -== — . (42) 

W 9+\n+x + 

The angular distribution of the intensity of the radiation field in Eq. (42) has the Lorentzian profile 
with the halfwidth 

A# = \lm X+ \. (43) 
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(9 

Figure 3. (color online) Functions \s(9 + n)\ (solid line) and s\(9 - n)\ (dashed line) for^+ = 0.1(1 - i). The 
width of the narrow peak at 8 = -? in the plot of \s(9 + n)\ is of the order of lmx+. For values of^f+ relevant 
to the THz range of frequencies, the peak is much narrower and higher. 

Similar calculations were done by V. Zon in [9, 10]. Note, however, that the range of applica- 
bility of Eqs. (40)-(42) is limited to very large distances, such as 

k r » \lm X+ r 2 - \s\. (44) 

For smaller distances, the factor s(8 - p) in Eqs. (39) cannot be considered smooth as compared 
with the exponent exp(-^ 2 ) since the former has a pole at p = 6 - |tt +^+, which approaches the 
saddle point p = —n at 6 = —\n. The absolute magnitudes of the functions s(G + n) are plotted in 
Fig. 3. More accurate calculation can be performed using the following approximate expression 
for the kernel function 

s{9 - p) * So (9 -p)= - 3 1 . (45) 

It is derived by expanding s(6 - p) in the Laurent series about the pole p but it provides a good 
approximation for s(6 - p) in the entire range of integration as Fig. 4 proves. Putting sq(6 - p) in 
Eq. (39) yields 

Bi(r, 6) = ±-.[s{e- p) e<^ dp = ~ f°° ^5 dt, (46) 

gi 
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Figure 4. (color online) Comparison of s(9 + n) (solid line) and Laurent expansion sq(9 + n) (dashed line); 
the real part is blue, the imaginary is black; x+ - 0.1(1 - i). 

where 

y/k r 1 + i ( n \ r — 

zo = j— (po + n) = — \6 + - +x+ 1 ihr. 

A method of calculating the integral 



W(z) 



2m J_ DO z - t 



in Eq. (46) is elaborated in the theory of plasma waves [22]. The result is expressed in terms of 
the imaginary error function 

erf 



rfi(z) = f exp (x 2 ) dx 
Jo 



and has two branches: 

W ± (z) = -- exp (-z 2 ) [+1 + i erfi(z)] . (48) 

The first branch W+(z) is originally computed in the assumption that the imaginary part of z is 
positive, Im(z) > 0, and then analytically continued to the lower half of the complex plane z, where 
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Im(z) < 0, using Eq. (48). The second branch W_(z), on the contrary, is originally computed for 
Im(z) < and then analytically continued to Im(z) > 0. In terms of the contour integration in 
Eq. (47), the analytical continuation implies that the pole t = z never crosses the contour of the 
integration (which originally goes along the axis lm(t) = 0), and the contour is deformed when the 
pole t = z crosses the axis lm(t) = to bypass the pole. Since the pole can be bypassed either from 
above or from below, there appear two branches of the function W(z). They differ in the residue of 
the integrand at the pole t = z: 

, 1 £ exp(-? 2 ) 
W + (z) - W.(z) = -exp(- z 2 ) = — (h PV dt. 

2m J z-t 

\t-z\=S 

For large \z\ » 1, the functions W+{z) are evaluated as 

1 1 , 

W + * — — , and W_ w — — + exp(-z 2 ), 
2 ymz 2 ymz 

respectively, if Im(z) > and as 



1 . . _. 1 



W + w — — - exp(-z z ), and W. 



2 ^fniz ' 2 yfiiiz 

if Im(z) < 0. Note that the exponent exp(-z 2 ) becomes very big in sectors where \n < \ arg(z)| < 



The sign of Im(z ) is reversed as the observation point at the angle 6 crosses the boundary 
6 w -\n of the shadow region. Together with the discontinuous field B , which represents the 
contribution of Eq. (34) of the pole p in the approximation of geometrical optics, the diffracted 
field B Y forms a continuous field. This total continuous field is given by the integral of Eq. (47) in 
Eq. (46) computed for Im(z) < 0. Indeed, the pole p goes away from the interior of the contour g 
through the contour gi if the angle 6 takes a sufficiently large negative value (which corresponds to 
a deep shadow region) thus making Im(z) < 0. In that case, the integral over the contour g gives 
no contribution to the result of calculations so B\ stands for the total field. Taking the lower sign 
in Eq. (48), we obtain the final expression 

B(r,9) = -W-\^-(6+ n - + X+ ) VVje^ (49) 

for the radiation field in the proximity of 9 = -\n. In plasma physics, our choice of the branch 
W-(z) of the multivalued function W(z) is known as the Landau rule [22]. By analytic continuation, 
Eq. (49) represents the total field for any values of 9. It is valid for any r larger than the wavelength 
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Figure 5. (color online) An angular pattern of a surface wave diffracted by a pure gold wedge for A = 140 /im 
and r = 1 mm (blue), r = 10 mm (magenta), and r = 50 mm (dark yellow). The x and y axes go along the 
upper and side faces of the wedge as shown in Fig. 1. The distance to the curve is proportional to the 
intensity of scattered radiation at the given angle. 

in a free space, r » 2n/kQ, both in the near-field and wave zones. The near field was computed 
numerically in [9, 10] by evaluating the integral in Eq. (14). Eq. (49) gives the desired result with 
a better accuracy, especially in case where Im(jf + ) — » 0. 

Angular patterns of diffracted waves computed by formula in Eq. (49) for a pure gold substrate 
and a gold wedge coated with ZnS film are presented in Figs. 5 and 6, respectively, where the 
intensity \B\ 2 is depicted for several distances from the edge of the wedge for parameters relevant 
to ongoing experiments at Novosibirsk Free Electron Laser. Fig. 5 is typical for the near- field zone, 
where the angular distribution has many peaks. Fig. 6, on the contrary, shows angular distributions 
with a single peak characterized by the Lorentzian profile of Eq. (42). In the first case of a pure 
gold wedge the far wave zone begins after the distance r ~ 15 m from the wedge edge whereas for 
ZnS -coated gold the wave zone is located at r « 20 mm. 

In contrast to the Lorentzian profile Eq. (42), the total diffracted field in Eq. (49) is not sym- 
metric about the boundary of the shadow region « It extends into the illuminated region 
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Figure 6. (color online) Same as in Fig. 5 for a gold wedge coated with Q.15jim ZnS film: r - 1 mm (blue), 
r = 10 mm (magenta, scaled x4), r = 50 mm and (dark yellow, scaled xl6). 

9 > ~n much further than into the shadow region 9 < - jn. This assertion has found a convincing 
evidence in experimental data cited in [3]. Another important conclusion is that the radiation in- 
tensity in the illuminated area does not copy the profile of the surface wave: instead of monotonic 
decreasing with the distance from the face of the wedge it rises at first. Again, this fact has been 
confirmed experimentally [3]. 

Intensity of diffracted waves vs. distance from the upper face is shown in Figs. 7 and 8, which 
illustrate that incident surface wave is mainly scattered into the upper hemisphere. 



VI. SCATTERED FIELDS NEAR WEDGE FACES 

Near the face surfaces of the wedge at 9 = the integrals B x and B 2 have commensurate 
magnitude and both should be computed. Moreover, the integrands in Eqs. (39) obey the equality 

s(9 + n) = s(9 - tt) (50) 

at 6 = ±\n as can be seen from Fig. 3. The approximation of Eq. (45) provides a good fit for the 
kernel s(9 - p) near the saddle point p = -n for almost entire interval -|tt < 9 < \n of the angles 
0. However it does not satisfy Eq. (50) at the ends and could be improved there. 

Near the upper face at 9 = |tt the pole p + = 9 -y+ - \n approaches the saddle point p = -n at 
9 = |7r, as mentioned at the end of Section IV. Expansion of the kernel s(9 - p) in Eq. (39a) into 
the Laurent series about the pole p+ = 9 - \n - x+ near the saddle point p = -n has the following 
form: 



s(9 -p)^ — — 1 + - n 

3V3V 9-ln- X+ -p) 
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Figure 7. (color online) Intensity of diffracted wave vs. vertical coordinate at different distances from the 
wedge edge: x = 1 mm (blue), x = 10 mm (magenta), x = 50 mm (dark yellow), and x = 150 mm (green). 
Pure gold wedge. 



The kernel 



s(9-p) * 



3V3 



(52) 



in Eq. (39b) is regular at 9 = |tt near the second saddle point p = n. Combining Eqs. (39), (51), 
and (52), we obtain 



Bi +B 2 



i r 

2ni J_ 



3V3 



1 + 



X+ 



2e ikor 



e— 2 ^ dt 
dt 



3 V3 



dp 

in/4 / 1 + i / 3 I r— 

X+ W. — [9--n-x + \^fhr 



y/2nk r 



■ (53) 



The second term in Eq. (53), containing is small and could be neglected. Eq. (49) gives a 
similar result at 9 = \n, which differs by a numerical coefficient of the order of unity. 

Similarly, the pole p = 9 + In + X- approaches the saddle point p = n at 9 — » -|tt. The 
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Figure 8. (color online) Intensity of diffracted wave vs. vertical coordinate at different distances from the 
wedge edge: x = 1 mm (blue), x = 10 mm (magenta, scaled x2), x = 50 mm (dark yellow, scaled x4), and , 
x = 150 mm (green, scaled x8). Gold wedge coated with a 0.75/mi ZnS film. 



1 (- 2 X _ 



expansion of the kernel of the integral in Eq. (39b) about p has the following form: 

1 - 

" V3 V 

and the kernel of Eq. (39a) is 

1 



0+ jn+x- -PJ 



s(9-p)<* -. 

V3 

Combining Eq. (39) with Eqs. (54) and (55) yields 

1 f 2 ( x- ) 
B l+ B 2 ^ — \ — 1 — 

_w/4 



(54) 



(55) 



dt 



2e ik ° r 



(56) 



Again, the second term with^_ is small and can be dropped. 
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Figure 9. (color online) Experimental configurations for the study of radiation induced by decoupled SPP 
using (a) a microbolometer focal plane array and (b) a Goley cell coupled to a lock-in amplifier SR830. 

Eq. (53) describes a wave propagating along the upper face of the wedge back from the wedge 
edge, and Eq. (56) is a wave propagating along the lower face. They can be thought of as reflected 
and refracted waves, respectively, although they are not surface waves. It is interesting that the 
amplitude of the refracted wave is 3 times larger than the amplitude of the reflected one. 

VII. COMPARISON WITH EXPERIMENT 

The experimental configurations are shown in Fig. 9. Monochromatic radiation of Novosibirsk 
Free Electron Laser (FEL) at a frequency of 2.3 THz (A = 140 /mi) entered the user station through 
a 16-m long beamline as a Gaussian beam / = Io exp (-2r 2 /w 2 ) with a waist of 9 mm. After 
passing a circular aperture with a diameter of 10 mm, the beam was focused with a cylindrical 
mirror into the input mouth of a plane waveguide, formed by the gold-covered facets of two glass 
prisms. At the output mouth of the waveguide, the radiation passed through the slit, transformed 
into a surface plasmon-polariton travelling along the metal-dielectric interface. Since a portion of 
the radiation could be emitted as a free wave, to separate the SPP and the bulk wave we made the 
input facet of the large prism tapered with an angle of 13 degrees to the upper sample plane. The 
latter was 17 cm long and 4 cm wide. The edge between these two facets was smoothed for the 
purpose of decreasing the SPP radiation loss. 

The large and small facets of glass slabs were covered with a 1 /im thick gold layer, which was 
considered within this problem as a bulk metal since the skin-depth for the gold is much smaller 
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than 1 /mi. The SPPs, which were launched with the help of the waveguide, travelled along the 
large facet of the samples. Bare gold and gold covered with ZnS layers 0.1 to 3 jjm thick were em- 
ployed in the experiments. To study SPP characteristics we applied two non-invasive techniques 
for detection of electromagnetic (EM) radiation in the free space behind the end facet of the sam- 
ples. An optical system consisting of a TPX lens with / = 50 mm [30] and a microbolometer focal 
plane array [31, 32] was used for imaging of the EM radiation wavefronts. Intensity of the EM 
radiation at different distances h was scanned along the v-axis using an opto-acustic Goley cell 
[30] with an input slit 0.2 mm thick directed along the z-axis. 

The experimental results were presented in brief in the conference proceedings [3] and will be 
published in detail elsewhere [33]. In this paper we present only the data relevant to the contents 
of the theory, namely we describe characteristics of the electromagnetic field (EMF) that arose 
when an SPP reached the end of the surface. The main features of the EMF were as follows (see 
Figs. 10 and 1 1). (a) The intensity of EMF, in contrast to the SPP intensity, reached its maximum 
at a distance of 1 to 2 mm above the sample surface, depending on the ZnS layer thickness, and 
only then sloped down, (b) Oscillations were observed on the descending part of the distribution, 
(c) The width of the distribution decreased with the distance and the angular distribution became 
narrower. All these features are in reasonable agreement with the theoretical predictions. Some 
discrepancies between the distribution shapes on the slope recorded with the MBFPA and the 
Goley cell may be caused, for example, by field disturbance by the metal diaphragm at the cell 
input window. 



VIII. CONCLUSIONS 

In this paper, we have considered diffraction of a surface wave by a rectangular wedge with 
impedance facets using the Sommerfeld-Malyuzhinets technique. We have derived Eq. (49) that 
uniformly describes the diffracted field both in the near-field and far-field zones. It was used in 
Section V for analysis of diffracted radiation at various distances from the wedge. Main conclu- 
sions from this formula are confirmed by available data from ongoing experiments at Novosibirsk 
Free Electron Laser [3, 33]. 

First, we have shown that the total diffracted field expressed by Eq. (49) is not symmetric about 
the boundary of the shadow region and that the surface wave is scattered mainly into the upper 
hemisphere. It extends into the illuminated region much further than into the shadow region as 
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Figure 10. Images of spatial distribution of radiation near the end of a bare gold sample (above) and a sample 
covered with a l-/im ZnS layer (below), recorded with the microbolometer focal plane array (configuration 
of Fig. 9,a). 

depicted in Figs. 7 and 8. 

Another important conclusion is that the radiation intensity in the illuminated area does not fol- 
low a surface wave profile and at first increases instead of monotonic decreasing with the distance 
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Figure 11. EM wave intensity distribution (in arbitrary units) vs. angle for three distances in the experi- 
mental configuration of Fig. 9, b for a bare gold surface: (a) Goley cell signal for three distances from the 
sample end; (b) simulations. 

from the face of the wedge, as one might expect for the surface wave, which exponentially decays 
with the distance from the metal-air interface. Again, this fact has been confirmed experimentally 
in [3, 33]. 

We have confirmed the conclusion made in [9, 10] that in the wave zone the angular distribution 
of the scattered wave has the Lorentzian form with a width determined by impedance. However, 
we have noted that the wave zone for a wedge with a small surface impedance begins at very large 
distances, specified by Eq. (44). 
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